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ABSTRACT 

We take advantage of the Hubble Ultra Deep Field (UDF) data to study the restframe optical and ultra 
violet (UV) morphologies of the novel population of Distant Red Galaxies (DRGs). Six galaxies with 
J — Ks> 2.3 are found to Kg = 21.5, five of which have photometric redshifts Zphot ^ 2 , corresponding 
to a surface density of 0.9 arcmin“^. The surface brightness distributions of the Zphot ^ 2 galaxies are 
better represented by exponential disks than R®/'^-laws. Two of the Zphot ^ 2 galaxies are extended, 
while three have compact morphologies. The restframe optical morphology of the Zphot ^ 2 galaxies is 
quite different from the restframe UV morphology: all the galaxies have red central components which 
dominate in the NICMOS iLieo-band images, and distinct off-center blue features which show up in 
(and often dominate) the ACS images. The mean measured effective radius of the Zphot ^ 2 galaxies 
is (re) = 1.9 ± 1.4 kpc, similar (within the errors) to the mean size of LBGs at similar redshifts. All 
the DRGs are resolved in the ACS images, while four are resolved in the NICMOS images. Two of the 
Zphot ^ 2 galaxies are bright X-ray sources and hence host AGN. One of these galaxies is resolved in 
the ACS and NICMOS images, which means the AGN does not dominate its restframe UV-optical SED, 
while the other is unresolved in the NICMOS images and hence could have an AGN dominated SED. 
The diverse restframe optical and UV morphological properties of DRGs derived here suggest that they 
have complex stellar populations, consisting of both evolved populations that dominate the mass and 
the restframe optical light, and younger populations, which show up as patches of star formation in the 
restframe UV light; in many ways resembling the properties of normal local galaxies. This interpretation 
is supported by fits to the broadband SEDs, which for all five Zphot ^ 2 are best represented by models 
with extended star formation histories and substantial amounts of dust. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution — galaxies: 
formation 


1. INTRODUCTION 

The Faint InfraRed ExtraGalactic Survey (FIRES, 
Labbe et al. 2003) discovered a highly clustered (Daddi 
et al. 2003) population of Distant {zphot ^ 2) Red Galax¬ 
ies (DRGs) by selecting galaxies with J — Kg >2.3 (Franx 
et al. 2003). At similar restframe optical luminosities, 
DRGs typically have older stellar ages (1-3 Gyr), larger 
stellar masses (^IO^^A/q), more dust (Ay=l-3 mag) and 
higher star formation rates (SFR^lOOM 0 i/r“®) than the 
complementary population of high redshift galaxies se¬ 
lected using the Lyman break technique (Steidel et al. 
1999; Shapley et al. 2001, 2003; Forster-Schreiber et al. 
2004; van Dokkum et al. 2004). The observed properties 
of LBGs suggests that they are on average younger than 
the more evolved DRGs. 

Morphological studies of Lyman Break Galaxies (LBGs) 
using the Hubble Space Telescope (HST) suggest that 
they are compact (re ~ 2 kpc) and have very similar rest¬ 
frame UV and optical morphologies (Giavalisco et al. 1996; 
Lowenthal et al. 1997; Dickinson 2000; Ferguson et al. 
2004; Bouwens et al. 2004). Here we present the first 


space-based study of the restframe optical morphologies 
of DRGs. We take advantage of the Hubble Ultra Deep 
Field (UDF) NICMOS and ACS imaging of a 2.4 x 2.4 
arcmin^ area of the sky in the Chandra Deep Field South 
(CDFS), to study restframe UV to optical morphologies 
of DRGs. The unprecedented depth® of the data allow for 
detailed morphological studies, even of very low surface 
brightness components. 

In the following we assume a standard cosmology with 
Q,m = 0.3, Da = 0.7, and h = 0.7 (Spergel et al. 2003). We 
refer to the HST filters F435W, F606W, F775W, F850LP, 
FllOW and F160W as H 435 , Vqqq, 1775 , zsbo, Jiio and 
Hiqo- Magnitudes are in the Vega system, unless other¬ 
wise noted. 


2. SELECTION 

In addition to the public UDF data (Beckwith et al. 
2004; Thompson et al. 2004) we use a photometric cat¬ 
alog consisting of ACS data (in the H 435 , Veoe, *775 
and zgso-bands) taken as part of the GOODS survey. 


^ Based on observations made with the NASA/ESA Hubble Space Telescope, which is operated by the Association of Universities for Research 
in Astronomy, Inc, under NASA contract NAS5-26555 

^ Department of Astronomy, Yale University, New Haven, CT 06520-8101, email: toft@astro.yale.edu, dokkum@astro.yale.edu 
^ Leiden Observatory, Postbus 9513, 2300 RA Leiden, Netherlands, email: franx@strw.leidenuniv.nl, mariska@strw.leidenuniv.nl 
Steward Observatory, University of Arizona, 933 Cherry Avenue, Tucson, AZ 85721, email: thompson@cosmos.as.arizona.edu 
^ University of California at Santa Cruz, 1156 High Street, Santa Cruz, CA 95064, email: gdi@ucolick.org, bouwens@ucolick.org 
® The lOo- AB limits in the filters F435W, F606W, F775W, F850LP, FllOW, F160W are 29.2, 30.0, 29.7, 28.7, 27.0, 27.0., (Beckwith et al. 
2004; Thompson et al. 2004) 
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Fig. 1.— Gallery of the six galaxies in the UDF with 
Ks < 21.5 and J — Ks > 2.3. The hlters shown are (from 
left to right): 5435 + Veoe, *775 + ^ 850 , ^iio, Hieo- The size of 
the stamps is 3'.'6 x 3'.'6. North is up, east is to the left. The 
images have been normalized such that a galaxy with constant 
f\ would appear the same in all images. The galaxies appear 
regular in the NICMOS images, and show a range of sizes. 

combined with deep VLT/ISAAC data (in the J, H and 
ATg-bands) taken as part of the ESO Imaging Survey (EIS) 
(van der Wei et al. 2004). In the following we refer to this 
catalog as the GOODS/EIS catalog. 

We use the GOODS/EIS catalog rather than the UDE 
catalog for selecting galaxies as it includes the same J and 
Ks filters that have been used previously (Franx et al. 
2003). The additional depth of the UDF Jno-band ob¬ 
servation (relative to the ISAAG J-band observations) 
is less important as we are concentrating on relatively 
bright galaxies. Galaxies were selected using the criteria: 
Ks < 21.5, J — Ks > 2.3. Six galaxies in the UDF match 
these criteria. A postage stamp gallery of the sample is 
shown in Fig. 1. 

3. PHOTOMETRY AND PHOTOMETRIC REDSHIFTS 

Photometric redshifts were derived both from 1"8 diam¬ 
eter aperture photometry derived from Point Spread Func¬ 
tion (PSF)-matched AGS/NIGMOS UDF images, and 
from 2" diameter aperture photometry available in the 
GOODS/EIS catalog. The GOODS/EIS dataset was 



found to provide stronger constraints on the redshifts than 
the UDF data set, due to the inclusion of the ATg-band 
and the fact that the ISAAG J-band is redder and nar¬ 
rower than the NICMOS Jno band. Hyperz (Bolzonella 
et al. 2000) was used to derive photometric redshifts. To 
limit the degrees of freedom in the fits, we used four em¬ 
pirical templates (E, Sbc, Scd and Im) to derive the best 
fitting photometric redshifts. In Tab 1, we list the best 
fitting photometric redshifts and other properties of the 
sample. The median photometric redshift of the sample is 
2.1 ±0.7, in good agreement with the median photometric 
redshift found for DRGs in the HDFS (Franx et al. 2003). 
One of the six galaxies in the sample has a photometric 
redshift significantly lower than 2. The surface density of 
Zphot ^ 2 DRGs to Ks = 21.5 is ~ 0.9 arcmin”^, consis¬ 
tent with ~ 1.0 arcmin“^ to Kg = 21.7 found for a much 
larger sample in the MS1054 field (Forster-Schreiber et al. 
2004). Following Franx et al. (2003) and Forster-Schreiber 
et al. (2004), we next fixed the redshifts and fitted three 
Bruzual & Chariot (2003) models (BC): a constant star 
formation model (CSF) with dust, a “Tau” model (ex¬ 
ponentially declining star formation with a characteristic 
timescale of t =1 Gyr) with dust, and a Single instan¬ 
taneous burst model without dust. In Fig 2. we show 
SEDs of the six galaxies in the sample, along with the 
best fitting redshifted BC models (full black curves). It is 
interesting to note that all six galaxies are better fitted by 
the “Tau” model with ongoing star formation and signifi¬ 
cant amounts of dust, than the single burst model with no 
dust (shown in dashed red curves). 


4. MORPHOLOGIES 

The DRGs in our sample show a range of morphologies: 
in the NIR, three are extended and three are compact. Of 
the galaxies with Zphot ^ 2, three are compact and two are 
extended. The observed NIR (restframe optical) morphol¬ 
ogy of the Zphot ^ 2 galaxies look quite different from the 
observed optical (restframe UV) morphology, highlighting 
the importance of deep NIR imaging of high redshift galax¬ 
ies. In the NIR, all the Zphot ^ 2 galaxies have red central 
components which dominate the NICMOS iJigo-band im¬ 
ages, while distinct blue off-center features show up in (and 
often dominate) the ACS images. We verified that these 
differences were not due to the different PSF and noise 
properties of the ACS and NICMOS images. 

The most pronounced differences between ACS and 
NICMOS are seen in two galaxies: UDF-DRGl, which 
has a luminous red bulge which dominates the light in the 
NIR bands and a more extended, patchy blue disk which 
dominates the light in the optical bands; and UDF-DRG6, 
which in the NIR bands is a single compact source, but 
in the optical bands is composed of four close compact 
sources of comparable brightness. From the optical bands 
alone these two galaxies could easily have been classified 
as highly irregular galaxies, or galaxies undergoing major 
merging, while their NIR morphologies reveal that they 
are more likely relatively relaxed galaxies with patches of 
ongoing star formation. 

The extraordinary depth of the UDF data allow 
us to test whether the blue “blobs” seen in UDF- 
DRG 1, 3, 5 and 6 are due to substructure in the 
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Fig. 2.— SEDs (/a versus A) of the six galaxies in the sample. 
Solid circles are photometry from the GOODS/EIS catalog 
(normalized by fx in ifs-band). The full black curves are the 
best fitting (SB, Tau or CSF) BC models and the red dashed 
curves are the best fitting SB models (see text for details). 
The Tau model (with some dust extinction) provides the best 
fit to all six galaxies. The CSE models provide almost as good 
fits, but require more dust. The formal errors on the redshift 
from HYPERZ are in the range ~0.1-0.2, as are error estimates 
from monte carlo simulations. 
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Fig. 3.— Comparison of red galaxies selected here from their 
J — Ks color (large solid symbols), from their 7775 — T/ieo 
color and z^hot (large open symbols, Chen & Marzke 2004), 
and from their 2350 — RAC) color (small filled (pink) 

symbols, Yan et al. 2004). The colors and magnitudes of the 
galaxies are from the COODS/EIS catalog. Four of the Chen 
& Marzke (2004) and four of the Yan et al. (2004) galaxies 
are too faint to be detected in our 72s-band data. The three 
samples only have two galaxies in common (UDF-DRC4 and 
6 ) 


galaxies, or chance alignment with unrelated sources at dif¬ 
ferent redshifts. To do this we placed small (3 pix) aper¬ 
tures on the positions of the individual blobs (indicated 
on the 77435 + heoe images in Fig. 1) in the PSF-matched 
images and derived their SEDs and photometric redshift. 
These are all consistent with the blobs being substructure 
at the redshift of the main galaxies, although we note that 
the photometry in neighboring blobs are not fully inde¬ 
pendent. 

5. SIZES 

To better quantify the differences and similarities of the 
galaxies and their appearance at different wavelengths we 
modeled their surface brightness distributions using gal- 
fit (Peng et al. 2002). We htted a Sersic profile (con¬ 
volved with the PSF, estimated from a calibration star): 
I(r) = /(0)ea;p[&„(r/re)^/"] which has been shown to 
be a good representation of the surface brightness of a 
range of galaxy types from ellipticals (n = 4) to disks 
(n = 1). Leaving the radial shape parameter n and the 
effective radius Vf. free in the fitting resulted in a mean 
value of n = 1.4 ± 0.9. Since there is some degeneracy 
between n and in the hts, we next fixed n to 1, to 
be able to compare the derived sizes of the galaxies with 
each other and with other high redshift galaxies. The re¬ 
sulting sizes (listed in Tab. 1) range from O'.'l to 0'.'5. 
At NIR wavelengths (in the TLigo-band), the mean mea¬ 
sured circularized radius^ of the hve Zphot ^ 2 galaxies is 
(re)NIR = 0"23 ± 0"16. At the redshifts of the galaxies 
this corresponds to a mean physical size of 1.9 ± 1.4 kpc. 
At optical wavelengths most of the sizes are slightly larger, 
but less well defined due to the patchy light distribution. 


All the DRGs are resolved in the ACS images. To deter¬ 
mine whether the “compact” galaxies (UDF-DRG 3, 4 and 
6 ) are resolved in the NICMOS data we repeated the size¬ 
fitting analysis on three stars (two stars in the field and 
a calibration star), to derive the mean size of “bonafide” 
unresolved objects: {ve)stars = O'.'OG ± 0"03. From this we 
conclude that UDF-DRG4 is marginally resolved, while 
UDF-DRG3 and 6 are consistent with being unresolved. 
Varying the fitting regions, and using different stars to 
model the PSF result in little variation in the derived r^- 

The derived sizes can be compared to those of 2 ~ 3 
LBGs. Ferguson et al. (2004) hnd a mean size (re) ~ 2.3 
kpc, and Bouwens et al. (2004) hnd (re) ~ 1.8 kpc. As the 
luminosity range of the ZpUot ^ 2 galaxies in our sample 
{Mv = —21.4 to —23.3) is comparable to that of the LBG 
samples we can compare the sizes directly (see Trujillo et 
al. 2004). The mean size for Zptot ^ 2 DRGs derived here 
is similar within the errors to the mean size of LBGs. 

A possible caveat of this comparison is the fact that we 
are comparing DRG sizes derived in the restframe optical, 
with LBG sizes derived in the restframe UV. However, the 
difference between restframe UV and optical sizes derived 
here for our small sample of DRGs (see Tab.l) are much 
smaller than the galaxy-galaxy variation in sizes. 


6. AGN CONTRIBUTION 

Two of the DRGs in our sample are detected in the 
IMs Chandra Deep Field South X-ray catalog (Giacconi 
et al. 2002) to fx > 10“^® ergs cm“^s“^. The [0.5-2.0] 
kev band X-ray fluxes of UDF-DRG 1 and 3 {fx = 73 and 
81 X 10“^^ergscm“^s“^ respectively) are comparable to 


^ defined as re,c = re\/(l -I- e), where is the effective radius encompassing half the flux, measured along the semi-major axis, and e = 1 — a/b 
is the ellipticity of light distribution 
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TABLE 1 

Properties of distant red galaxies in the UDF 


Id“ 

Kg 

J-Kg 

z*’ 

phot 

j>C 

' e,c 

UDF-DRCl /41 

19.38 

2.42 

2.2 

0'.'35/0" 40 

UDF-DRC2 /372 

20.69 

2.32 

1.2 

0'.T7/0'.T9 

UDF-DRC3 /290 

20.70 

2.66 

1.9 

0'.'09‘^/0'.'03 

UDF-DRC4 /746 

20.93 

3.26 

3.4 

0'.T4/0'.T2 

UDF-DRC5 /899 

19.81 

2.45 

1.8 

0'.'49/0'.'53 

UDF-DRC6 /1207 

21.00 

3.24 

3.4 

O'.'08'^/O'.'IO 


“ Id in this letter / Id in the public NICMOS UDF catalog (vers.2.0) 

Photometric redshift derived from the GOODS/EIS catalog 

Circularized effective radius derived in the iLieo/Veoe filters. The 
typical error is ~ 0”1 
“ Consistent with unresolved 

the fluxes of the two DRGs with AGNs in the MS 1054- 
03 field (Rubin et al. 2004). The corresponding X-ray 
luminosities are too bright to be supported by star forma¬ 
tion and must thus be of AGN origin (the X-ray flux, due 
to star formation, expected from DRGs fx — 10“^^ergs 
s“^cm“^ (Rubin et al. 2004) is below the detection limit). 
The fact that two of the Zphot ^ 2 galaxies certainly hosts 
AGN allow us to test the AGNs contribution to the rest- 
frame optical emission. UDF-DRGl is clearly resolved in 
the NIR bands, so its red restframe optical color is prob¬ 
ably due to evolved stars, rather than the AGN. UDF- 
DRG3 on the other hand is unresolved in the NIR bands, 
and could thus be dominated by the AGN, rather than 
stars, in its restframe optical emission. 

7. DISCUSSION 

Two other studies have recently identified red galaxies 
in the UDF. Chen & Marzke (2004) select 9 galaxies with 
I 775 — Hiqq > 2 and Zphot > 2.5, six of which are ar¬ 
gued to be best described by evolved templates with no 
dust, and three by dusty starbursts. Yan et al. (2004) 
combine the HST data with mid-infrared observations ob¬ 
tained with the Infra Red Array Camera (IRAC) onboard 
the Spitzer Space Telescope, to select 10 galaxies with 
/^(3.6/im)//^(2;85o) > 20, and argue that their SEDs are 
well represented by a two component model consisting of 
an old evolved population and a secondary young popula¬ 
tion. In Fig. 3 we plot J — Kg vs Kg for the three samples 
of red galaxies in the UDF. 

The two reddest DRGs in our sample make it into both 
the Chen & Marzke (2004) sample (Object 09024 and 
01927 in their Table 1) and the Yan et al. (2004) sam¬ 


ple (object 8 and 9 in their Table 1). Of the remaining 
seven/eigth galaxies in the Chen & Marzke (2004)/Yan et 
al. (2004)) samples, four/four are too faint to be detected 
in our data, while three/four are either too blue or too 
faint to make our selection criterion. We conclude that 
while the different selection methods show some overlap, 
they are not redundant. 

We note, that while the photo-z’s derived for the two 
galaxies present in all three samples agree reasonably well 
(±0.2), the model parameters are uncertain: we find both 
galaxies to be best fit by a IGyr Tau model with dust; 
Chen & Marzke (2004) find UDF-DRG4 to be best fit by 
a 0.3Gyr Tau model without dust, and UDF-DRG6 by a 
dusty CSF model; and Yan et al. (2004) fits a two com¬ 
ponent model without dust to both galaxies. This demon¬ 
strates that while photometric redshifts derived from fits 
to broadband SEDs are relatively robust, parameters of 
the fitted models are in general less well determined and 
should be interpretted with caution. 

The most characteristic property of the restframe opti¬ 
cal and UV morphologies of the DRGs studied here is their 
diversity. Some are extended, some are compact, some 
have very different restframe optical and UV morphologies, 
while some look similar. The conclusion that can be drawn 
is that most of the DRGs in the present sample have com¬ 
posite stellar populations, with central populations of red 
evolved stars dominating the restframe optical morpholo¬ 
gies and younger more patchy distributed blue stars dom¬ 
inating the restframe UV morphologies, while LBGs have 
similar restframe optical and UV morphologies (Dickinson 
2000), and thus may be uniformly younger systems with 
star forming regions so luminous that they dominate even 
the restframe optical light. These results are in agreement 
with those derived from spectra and broadband SEDs of 
DRGs (Forster-Schreiber et al. 2004) and LBGs (Pettini 
et al. 2001; Shapley et al. 2003). There might, however, 
be some overlap in the properties of restframe optical and 
UV selected galaxies as a small fraction of restrame UV 
selected LBGs at z ~ 3 and Bm/Bx galaxies at z ~ 2 have 
recently been shown to have red optical colors and ages 
similar to the DRGs (Steidel et al. 2004; Shapley et al. 
2004). 

Support from the Danish Natural Science Research 
Council, and from NASA through HST grant GO-09803- 
10-A and LTSA NNG04GE12G is gratefully acknowl¬ 
edged. 
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